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ABSTRACT 

In this work, we have stndied the observational constraints on the Cardassian 
Model for the dark energy. We have compared the model with existing Snpernova 
data. The dependence of the locations of the Cosmic Microwave Backgronnd 
Radiation (CMBR) peaks on the parameters of the model have also been stndied. 

We hnd, in particnlar, that observational data arising from Archeops for the 
location of the hrst peak, BOOMERANG for the location of the third peak, 
together with the Snpernova data, constrain signihcantly the parameter space. 

Subject headings: Dark Energy, Cardassian Expansion, CMB, SNIa 


1. Introduction 

It is remarkable that number of current observations indicate that we are living in 
a spatially flat, low matter density universe which is currently undergoing an accelerating 
expansion (Bernadis et al. 2000; Hanany et ah 2000; Balbi et al. 2000; Perlmutter et ah 1997, 
1998; Garnavich et al. 1998; Riess et al. 1998). The most simple explanation of this current 
cosmological state of the universe requires two dark components: one is in the form of non- 
relativistic dust (“dark matter”) with vanishing pressure contributing one-third of the total 
energy density of the universe and clustering gravitationally at small scales while the second 
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one is a smoothly distributed component having large negative pressure (“dark energy”) 
and contributing around two-third of the total energy density of the universe. Although the 
simplest candidate for this dark energy is the vacuum energy or the cosmological constant 
(A), alternative scenarios where the acceleration is driven by dynamical scalar field both 
minimally (Caldwell et al 1998; Peebles et ah 1988; Ferreira et ah 1987; Copeland et ah 
1988; Steinhardt et al. 1999; Zlatev et al. 1999; Wetterich 1988; Ratra et al. 1988; Barreiro 
et al. 2001; Sahni et al. 2000; Sen et al. 2002; Bento et al. 2001) and non-minimally (Bertolo 
et al. 1999; Bertolami et al. 2000; Uzan 1999; Amendola 2000; Gasperini 2001a; Sen et 
al. 2001; Sen and Sen 2001a,b) coupled with gravity called “quintessence” have been widely 
investigated in recent years 

As none of the two components (dark matter and dark energy) has laboratory evidence 
both directly or indirectly, one has to invoke untested physics twice to explain the current 
observations. That is why people in recent times have proposed interesting scenarios where 
one describes both dark matter and dark energy in a unified way through a single fluid com¬ 
ponent in the Einstein’s equation. Chaplygin gas model is one such interesting possibility 
which has attracted lot of attentions in recent times (Kamenshchik et al. 2001; Bento et al. 
2002a,b; Bilic et al. 2002). Padmanabhan and Roy Choudhury have also proposed an inter¬ 
esting unified description based on a rolling tachyon arising in string theory (Padmanabhan 
et al. 2002). 

Recently Freese & Lewis(Freese et al. 2002) proposed another interesting alternative 
to quintessence scenario where the recent acceleration of the flat universe is driven solely 
by the matter, instead of using any cosmological constant or vacuum energy term. Since 
pure matter or radiation cannot alone take into account the recent acceleration in the flat 
universe, this goal is accomplished by modifying the Friedman equation with an empirical 
additional term named Cardassian term. 

H^ = Ap + Bp^, ( 1 ) 

where A = and B and n are constants and are the parameters of the model. Here the 
energy density (p) contains only matter (p^) and radiation (p^), i.e, p = Pm + Pr- Since at 
present pm» Pr, P can be considered consisting of pm only. The new term, dominates only 
recently at redshift ~ 1. To provide the required acceleration of the universe as the outcome 
of the dominance of this term, n should be < 2/3. 


^See also (Banerjee et al. 2001a; Banerjee et al 2001b; Sen and Seshadri 2000; Amendola 1999; Chiba 
1999; Gasperini 2001b; Riazuelo et al. 2001; Tocchini-Valentini et al. 2001) 
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There are several interpretations for the origin of this new “cardassian term” appearing 
in the Einstein’s equation (1). As described in (Freese et ah 2002) and also in (Chung et 
ah 2000), this term can appear as a consequence of embedding our observable universe as 
a 3 + 1 dimensional brane in extra dimension. Although recently it has been argued (Cline 
et al. 2002) that cardassian model based on this higher dimensional interpretation violates 
the weak energy condition for the bulk stress energy for n < 2/3 which is necessary for 
accelerating universe in late times. This extra term may also arise due to the matter self¬ 
interactions that contributes a negative pressure, through a long-range conhning force which 
may be of gravitational origin or may be a hfth force (Gondolo et al. 2002). Also denoting 
the second term as one can recast equation (1) as, 

H^ = A{p + p,), (2) 


where one can consider px as the dark energy with equation of state 

Wx = {n-l) + - - 

3 “1“ PmO^ 

where p^o and pmo are the present energy density for radiation and matter and we have as¬ 
sumed that present scale factor oq = 1. At late times, when pr « Pm, one can approximate 
the above equation of state as 

u>“' = (n-l) + /^(l+2) (4) 

^ PmO 

which is very slowly varying function with redshift. This is the conventional first order 
expansion to the equation of state for the dark energy and has been widely used in litera- 
ture(Weller et al. 2001; Astier 2001; Linder 2002). Also as pro « Pmo, at late times, 
is almost constant and it is identical to a dark energy component with a constant equation 
of state(Avelino 2002). But in early times, as one can not ignore the radiation component, 
one has to take the general equation of state Wx which is not constant. This is important 
when one considers the constraints on the model from the CMBR observations. 


In this work, we shall consider the observational constraints on the different parameters 
of this cardassian model. We shall consider the constraints arising from the positions of the 
peaks of the CMBR as well as those arising from Supernova observations. 


2. Cardassian Model 

One can cast equation (1) in the following way 

H^ = A p[l + (5) 

pear 
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where pcar is the energy density at which the two terms are equal. Once the energy density 
p drops below p^ar the universe starts accelerating, p^ar is given by 


pear 



Pmo(l + ^car)^{l + -(1 + ^car ) }■ 

^ ‘mO 


( 6 ) 


where Zcar is the redshift at which the second term, in equation (1), starts dominating over 
the hrst term. The model has two main parameters B (or pcar or Zcar) and n. 


To fulhll the requirement of the CMBR observation of a flat universe one can modify 
the critical energy density pc so that the matter can be sufficient to provide a flat geometry. 
Evaluating equation (1) today, 


B 

Hq = A {pmO + Pro){l + (PmO + Pro)^ 


( 7 ) 


In the new picture, Oq is dehned as so that matter alone makes the geometry flat. 


pc 

Here the expression for the critical energy density pc has been changed from its usual one 

72 , 

), as. 

Pc = pea X E(n, Zcar) (8) 


Pca{= ^), as. 


where 


B 


F{n,Zcar) = [1 + ^Pmo'(l + ^ 

i imO 

= |l + (l + j,„pi-”){l + ^(l + z,„)}>-”(l + ^)'-]-> (9) 

i LrnO ^ ^mO 

Qfno and are two parameters dehned as and hlj-o = ^ respectively. 


So the new critical density pc is expressed as a function of two parameters n and B or 
Zcar- This is similar to the expression given in (Freese et al. 2002). As mentioned there, 
a point to note here is that the new critical density pc is now the fraction of the original 
critical density pca which has a standard value 1.88 x 10“^® /iq gm/crrA. Hence the new 
critical density can be much lower than the standard estimate. And also keeping in mind 
that we consider a hat geometry (Hq = 1), we have today’s energy density p^o + Pro = Pc he, 

+ hlrO = = F (10) 

Pca 

In hgure 1, we present the diherent combinations of n and Zcar for certain values of F.However, 
gravitational cluster(Calberg et al. 1996) and other data suggest (Turner 2001) the total mat¬ 
ter density to be 30% of the usual critical density i.e, pc = .3 pca- This sets a preferred value 
.3 for F. 
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Now substituting the evolution of matter and radiation, we write equation (1) as the 
following, 

= A [pm^a ^(1 + T ® 

i imO ^ ^ ‘mO 

From equation (9) it is very straight forward to express B in terms of and hlmo 


|Ao' = (- 


^rO ^ ^rO \ —r, 


n ' a ' ■ 

a LffiQ i ijTiO 

Substituting this expression in equation (11), one can hnally recast equation (1) in the 
following fashion 


= QmoHQtt 


/ . ^rO \ I 

(a + —) + a 


—4n+4 


n 


mO 




mO 


a 


mO 


1 + ^ / I 

^'mO / J 


(13) 


In equation (13) is expressed in terms of the two model parameters n and hlmo which, in 
turn, is related to the other form of this parameter (F or Zcar) through equation (9) and (10). 
We are now in a position to constrain these two parameters with different observations. Once 
we constrain these two parameters from the observation we can evaluate the corresponding 
value of F or Zcar from eqn (10) and hgure 1. There is an important point to note is that 
the value n = 0 corresponds to a ACDM model (O^ at present is negligible). This will be 
crucial when we shall talk later about the allowed region of the parameter space. 


The prior assumptions in our subsequent calculations are as follows: scale factor at 
present oq = 1, scale factor at last scattering ais = 1100“^, h = 0.65, density parameter for 
radiation and baryons at present 11^0 = 9-89 x 10“^, hlbo = 0.05, and spectral index for the 
initial energy density perturbations, n = 1. 


3. Fitting with Supernova data 

First we check the consistency of the model with SNIa observations. The data from 
different supernova observation is related to a quantity called luminosity distance di dehned 
by 

di = {l + z)ri (14) 

for a source at r = ri at t = ti. But basically the luminosity distance in logarithmic units is 
what is observed by the astronomers. 


mB{z) = M + 51og,|,(D() 


(15) 
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where M. = M — 5 logj^Q Hq + 25 and Di = Hodi is the dimensionless luminosity distance. To 
measure Ai, one can show that for nearby sources (in low redshift limit) the above equation 
can be approximated as 

msiz) = M + 5 logio(2:). (16) 

The low redshift supernovae measurements can be used to calculate A4. Using equation 
(14), we estimate the magnitudes of the supernovae at different redshifts from 


tub- M = 51ogio{ 


where 




(1 + ;^) r dz' 

Jo {l + ZyX{z) 


} 


(17) 


a 


+ ^) + (i + ^) 


l + z fl. 


4n—4 


mO 




mO 


n 


mO 


1 I ^rO 

1 + 2 UmO 

1 _|_ ^rO 
^mO 


Using the above the relation to estimate ms at different redshifts and the observed 
values of the effective magnitude and the same standard errors and cr e// for a 

given redshift as listed by Perlmutter et al from SCP(Perlmutter et ah 1997), we compute 
as 


X 


2 


E 


- ruBizi)]"^ 
(^z,i + 


(18) 


We consider the data set of 54 supernovae comprising of 38 high redshift supernovae from 
Supernova Cosmology Project together with 16 low redshift supernovae from the Calan- 
Tololo project, as used by (Perlmutter et ah 1997) in their primary £t C (for details of 
excluded data points see (Perlmutter et al. 1997)). In figure 2 we present the permitted 
parameter space by the supernova constrains at different confidence level. We observe that 
the model best fits the current supernova data at 80% confidence level. 


4. Constraints from CMBR 

Our second tool for constraining the parameters is CMBR anisotropy spectrum. The 
CMBR peaks arise from oscillation of the primeval plasma just before the universe becomes 
transparent. The oscillation of the tightly bound photon-baryon fluid is a result of the 
balance between the gravitational interaction and photon pressure and this oscillations gives 
rise to the peaks and troughs in the temperature anisotropic spectrum. In an ideal photon- 
baryon fluid model, there is an analytic relation for the location of the m-th peak(Hu et al. 
2001 ): 




(19) 
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where Ia is the acoustic scale which depends on both pre and post recombination physics 
and also on the geometry of the universe. This has an analytical expression given by ^ 
where D is the angular diameter distance to the last scattering and sis is the sound horizon 
at the last scattering. In terms of the conformal time r, I a is given by (Doran et al. 2001, 
2002a), 


/a = TT 


^0 '^Is 


( 20 ) 


where tq and tis are the conformal time today and at last scattering and Cg is the average 
sound speed before last scattering. Cg is a constant for a particular —. We take it as 0.52 as 
others(Doran et al. 2001). 


Now, to hnd Ia, we write equation (13) in terms of conformal time 


,da, 

'dr' 


_ Q tT2 


/ I ^rO \ I 

(a + ——) + a 


—4n+4 


n 


mO 


1 




a 


mO 


(X -\- 


^mO 


1 _|_ ^rO 


( 21 ) 


wherefrom it is quite easy to hnd 




and 


Tls = 


To = 


dr = 


dr = 




mO-^O -^0 

1 r 


da 

X{a) 


where X{a) = ^{a + 


da 

alSnJo W) 


( 22 ) 


(23) 


a+ 


^mO 


1+7 


Substituting the above expression in equation (20), we have the analytical expression 
for IA in case of this model 

fl da 

Jo 


Ia = 


TT 


/c 


/o X{a) 
^is da 


- 1 


(24) 


^ [Jo X(a) 

where aig = 1100“^. Thus we can hnd the positions of the peaks in the CMBR spectrum 
from equation (19). I a and consequently the positions of the CMBR peaks here depends only 
n and Dmo- So once we calculate the positions of the peaks we can constrain the parameters 
n and Dmo by comparing the results from diherent observations. 

Now the simple relation in (19) is modihed by driving and dissipative ehects which 
introduces a phase shift to the oscillation such that(Hu et al. 2001) 


— ^Aij^ 0m) 


(25) 



The phase shift of the peaks (pm is predominantly determined by the pre-recombination 
physics and is independent of the geometry of the universe. It depends on parameters like, 
n, and the ratio of the energy density of radiation to matter at last scattering. For 
n = 1 and = 0.02 and rig = ~ 1100 at last scattering, the phase shift for the 

hrst peak(Hu et al. 2001) 

01 = 0.267(^)'b (26) 

Now substituting for ris a.t z = 1100, we have 


= 0.267( 


1100^1^o^ 

0.3flmo 


(27) 


Hence the position of the hrst peak /i is 


h — ^u(l — 0i) 


(28) 


Using equation (22)-(24) and (27) and (28) one can calculate li as a function of n and flmo. 
The observational bounds on li as predicted by BOOMERANG(Bernadis et al. 2002) and 
more recently by Archeops (Benoit et al. 2002) are li = 221±14 and li = 220±6 respectively. 
As the bound coming from the Archeops data is more stringent, we shall take this bound 
for hrst peak to constrain our parameters. 

The relative shift of the second peak is a very sensitive quantity and depends on many 
parameters. Hence it is very difficult to derive any constraint from the second peak. So we 
disregard the second one. As far the third peak is concerned Doran et al(Doran et al. 2002b) 
have shown it to be insensitive to diherent cosmological parameters. They estimated 03 to 
be 0.341. Hence 

/3 = /^(3 - 03) = /^(3 - 0.341) = 2.659 U (29) 

With this expression one can also calculate I 3 for diherent values of n and Dmo- The observa¬ 
tional bound on as suggested by BOOMERANG : /s = 845^25- hgure 3 we have shown 
the constraints on the parameter space that are obtained from the observational bounds on 
the location of the hrst (dashed contour) and third(full contour) GMBR peaks. Hence, from 
the GMBR point of view the allowed region of the model parameters lies in the intersection 
between these two contours. 


5. Conclusion 

In conclusion, we have shown that the locations of the GMBR peaks, as determined 
via Archeops and BOOMERANG, as well as the present Supernova la data, constrain a 
sizable portion of the parameter space of the Gardassian model. We observe that the model 



best fit with the Supernova data at 80% conhdence level. This, together with the allowed 
region from CMBR data restricts the parameters as 0.31 ^ n < 0.44 and 0.13 < hlmo ^ 0.23 
(See hgure 4). This clearly does not include the n = 0 which is the corresponding ACDM 
case in this model. This is the most interesting result of this investigation. Also the joint 
analysis indicates a lower value of flmo- This is consistent with that predicted by Zhu 
et.al (Zhu 2002) by investigating the constraints on the cardassian model from the recent 
measurements of the angular size of high-z compact radio sources. Also our bound on 
is consistent with that predicted by Melchiorri et.al (Melchiorri 2002) from a combined 
CMB+HST+SNIa+2dF analysis. Clearly with future high precision measurements of the 
MAP and PLANCK mission, we expect that the positions of the CMBR peaks will be 
determined with higher accuracy. This, together with the upcoming data from future SNAP 
mission will further constrain the parameter space of this model. 
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